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Abstract 



According to general relativity, the gravitomagnetic Lense-Thirring 
force of Mars would secularly shift the orbital plane of the Mars Global 
Surveyor (MGS) spacecraft by an amount of 1.5 m, on average, in the 
cross-track direction over 5 years. The determined cross-track post- 
fit residuals of MGS, built up by neglecting just the gravitomagnetic 
force in the dynamical force models and without fitting any empiri- 
cal cross-track acceleration which could remove the relativistic signal, 
amount to 1.6 m, on average, over a 5-years time interval spanning 
from 10 February 2000 to 14 January 2005. The discrepancy with the 
predictions of general relativity is 6%. 
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1 Introduction 

In this communication we propose to explain certain measured features of 
the path of the Martian orbiter Mars Global Surveyor (MGS) (Konopliv et 
al. 2006), whose orbital parameters are quoted in Tabled with the action of 
the general relativistic gravitomagnetic force on the spacecraft. If confirmed 
by subsequent, dedicated analyses it would be, in general, the first test of 
general relativity in the gravitational field of a Solar System body other 
than the Sun and the Earth and, in particular, the first independent test 
of the Lense-Thirring effect with respect to those conducted so far with the 
LAGEOS satellites in the Earth's gravitational field (lorio and Morea 2004, 
Ciufohni and Pavhs 2004, lorio 2006). 



1 



2 The Lense-Thirring effect 



The Lense-Thirring effect consists of tiny secular precessions of the longitude 
of the ascending node Q and the argument of pericentre w of a test particle 
freely orbiting a central rotating body. They are (Lense and Thirring 1918; 
Barker and O'Connell 1974; Soffel 1989; Ashby and Allison 1993; lorio 2001) 
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where G is the Newtonian gravitational constant, S is the spin angular 
momentum of the central body which acts as source of the gravitational 
field, c is the speed of light in vacuum, a, e and i are the semimajor axis, 
the eccentricity and the inclination of the orbital plane to the central body's 
equator, respectively. In general, the gravitmagnetic force affects the in- 
plane along-track and the out-of-plane cross-track parts of the orbit of a the 
test particle. Indeed, the shifts along the radial R, along-track T and cross- 
track N directions can be, in general, expressed in terms of the Keplerian 
orbital elements as (Christodoulidis et al. 1988) 
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where A4 is the mean anomaly. The Lense-Thirring shifts are, thus 

r Ai^LT = 0, 



(2) 



J ATlt = a-Jl + % (Awlt + coszAJ^lt) , 



(3) 



AA^LT = 0^/1 + V sin zAJ^LT- 



For a near polar satellite like MGS only the out-of-plane shift AiVLx is 
present. 



2 



Table 1: Orbital parameters of MGS and its Lense-Thirring precession, a 
is the semimajor axis, i is the inclination, and e is the eccentricity. J^lt is 
the gravitomagnetic node precession and (AA^lt) is the average out-of-plane 
gravitomagnetic shift over five years. 

a i e ^LT (AA^lt) 

3792.42 (km) 92.86 (deg) 0.0085 33.8 (mas yr"^) 1.5 (m) 



3 The interpretation of the MGS data 

In (Konopliv et al. 2006) six years of MGS Doppler and range tracking data 
and three years of Mars Odyssey Doppler and range tracking data were 
analyzed in order to obtain information about several features of the Mars 
gravity field summarized in the global solution MGS95J. 

As a by-product of such an analysis, also the orbit of MGS was deter- 
mined with great accuracy (Konopliv et al. 2006). Here we are interested in 
particular in the out-of-plane portion of its orbit. The root-sum-square out- 
of-plane residuals over a five-years time interval spanning from 10 February 
2000 to 14 January 2005 are shown in Figure ^ They have been determined 
in a Mars-centered coordinate system defined by the inertial International 
Celestial Reference Frame (ICRF), which is nearly equivalent to the Earth's 
mean equator at the epoch J2000, by using all the previously estimated pa- 
rameters of the global solution. Due to improved modelling (orientation, 
gravity, angular momentum wheel de-saturations, and atmospheric drag), 
the average of such residuals amounts to 1.6 m 

This result can be well explained with the action of the gravitomagnetic 
field of Mars on the orbit of MGS. Indeed, by using for the spin of Mars 
the value^ S = 1.9 x 10^^ kg m^ s~^, the Lense-Thirring precession of the 
node of MGS amounts to^ Clu^ = 33.8 milliarcseconds per year (mas yr~^), 
which yields a shift in the out-of-plane direction AA'lt = asinzAr^LT = 0.6 
m yr~^. This means that over five years the orbital plane of MGS was 

^It comes from S = aMB?{2-K/T), where the polar moment is a = C/MR^ = 0.3654 
(Konophv et al 2006) and the sidereal rotation period is T = 1.02595675 d (Seidelmann 
1992). 

■^The formulas of eq. (0 hold in a frame whose reference plane is the central body's 
equator. This is just the case because the Mars equator is tilted only by 1.7 deg to the 
Earth's equator, which is the reference plane of the frame used in the analysis of (Konopliv 
et al. 2006). Indeed, the obliquity of Mars and the Earth amount to 25.2 deg and 23.5 
deg respectively. 
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Figure 1: Orbit overlap differences (10 February 2000-14 January 2005) 
for the out-of-plane direction (Data provided by A. Konopliv, 2006). The 
average amounts to 1.6 m: the average Lense-Thirring out-of-plane shift, 
over the same time span, is 1.5 m. 
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shifted by 3 m by the gravitomagnetic field of Mars, a quantity which falls 
within the reached orbit accuracy. The averaged Lense-Thirring shift over 
five years is thus 1.5 m. 

To corroborate this interpretation of the data, it is important to note that 
in the suite of dynamical force models used to reduce the gathered observa- 
tions general relativity is, in fact, present, but only with the gravitoelectric 
Schwarzschild part: the gravitomagnetic Lense-Thirring part is not included, 
so that it is fully accounted for by the showed residuals. Moreover, no em- 
pirical out-of-plane accelerations were fitted in the data reduction process; 
this was done only for the along-track component of MGS orbit to account 
for other non-gravitational forces like thermal, solar radiation pressure and 
drag mismodelling. Thus, it is unlikely that the out-of-plane Lense-Thirring 
signal was removed from the data in the least-square procedure. Moreover, 
if the relativistic signature was removed so that the determined out-of-plane 
residuals were only (or mainly) due to other causes like mismodelling or un- 
modelling in the non-gravitational forces, it is difficult to understand why the 
along-track residuals (middle panel of Figure 3 of (Konopliv et al. 2006)) 
have almost the same magnitude, since the along-track component of the 
MGS orbit is much more affected by the non-gravitational accelerations 
(e.g. the atmospheric drag) than the out-of-plane one. 

4 Conclusions 

In this communication we showed that the general relativistic gravitomag- 
netic Lense-Thirring effect can accommodate 94% of the 1.6 m average error 
in the determined out-of-plane orbit residuals of the Martian or biter MGS. 
The collection and the processing of more data from MGS, Odyssey and 
the recently launched MRO spacecraft will be of the utmost importance in 
corroborating such test. For example, the Doppler tracking data from the 
Mars Exploration Rover (MER) to Odyssey will improve the determina- 
tion of just the out-of-plane portion of the Odyssey orbit. An independent 
check which could, and should, be done is to produce a particular solution 
in which it is also introduced a parameter which explicitly accounts for the 
Lense-Thirring eff'ect and solve for it. 
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